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I. INTRODUCTION 

The solution of the Schrodinger equation with position-dependent mass for any spherically 
symmetric potential has attracted attention over the past years [1-17]. The motivation in this 
direction arises from considerable applications in the different fields of the material science 
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and condensed matter physics. For instance, such apphcations in the case of the bound 
states in quantum system [4], the nonrelativistic Green's function for quantum systems with 
the position-dependent mass [5], the Dirac equation with position-dependent mass in the 
Coulomb field [5], electronic properties of semiconductors [10], ^He cluster [11], quantum 
dots [12], semiconductor heterostructures [13,14], quantum liquids [15], the dependence of 
energy gap on magnetic field in semiconductor nano-scale quantum rings [16], the solid 
state problems with the Dirac equation [17]. Almost all of those works mentioned above were 
focused on obtaining the energy eigenvalues and the potential function for the given quantum 
system with the position-dependent mass. The wave functions were either obtained by the 
solutions to the Schrodinger equation with the constant mass, or a few lower excited states 
were obtained by acting of the creation operator on the ground state. The effective potentials 
are the sum of the real potential form and the modification terms emerged from the location 
dependence of the effective mass [2] . Taking into consideration the position-dependent mass, 
the aim of this work is to carry out the analytical solutions of the A^— dimensional Schrodinger 
equation with position-dependent mass for a general class of static quarkonium potentials 
by the series expansion method used in [18,19]. Additionally, we investigate the linear-plus- 
Coulomb (Cornell) potential case [20,21]. 

The contents of this paper is as follows. In Section II, we present the A^ dimensional 
Schrodinger equation with position-dependent mass for any spherically symmetric potential. 
In Section III the analytical bound-state series solutions of a general class of quarkonium 
potentials. In Section IV, we study the analytical series solutions of the Cornell potential 
for a particle of an exponentially decaying mass, m = rriQe^^^, A > 0, case and then give an 
example for the analytical calculations of the Coulomb's wave function. Finally, in Section 
V we give our concluding remarks. 



II. THE AT-DIMENSIONAL POSITION-DEPENDENT MASS SCHRODINGER 

EQUATION 

The wave Schrodinger equation with position-dependent mass for a spherically symmetric 
potential V{r) in A^-dimensional space {in h = 1 units): 

Vjv^ Vjv ^(r) + 2[E- Vir)] ^{r) = 0, (1) 

where the wave function is defined by [20,21,22,23] 

^(r)=e-(^-i)/2i?„,Kr)n„(x). (2) 



and m = m{r). We substitute Vat^Vat V^(r) = (Vat^J • (Vn i'{r)) + ;;^V^ V^(r) into 
Eq.(l) and obtain the A^— dimensional position-dependent mass radial Schrodinger equation 

[ dr^ -m y 2r dry 4r^ J 

where k = N + 21 and m'{r) = dm/ dr. It can be clearly seen that for m' = case, the above 
equation reduces to the well known equation with constant mass used in Refs. [20,21,22,23]. 
In the present work, we are concerned in bound states, i.e., E < 0. On the other hand, one 
should be careful about the behavior of the wave function R{r) near r = and r — > oo. It 
may be mentioned that R{r) behaves like r*^^^^)/^ near r = and it should be normalizable. 
We choose the wave function [18,19,24,25] 

Rn,i{r) = r('^-^)/2e-^'-M(r), b = ^J-2moE, (4) 

where mo is the initial value of mass. It should be noted that N and / enter into expression 
(3) in the form of the combination k = N + 21. Consequently, the solutions for a particular 
central potential V{r) are the same as long as k remains unaltered. Thus, the s-wave 
eigensolution (Rno) and eigenvalues (E) in four-dimensional space are identical to the p- 
wave two-dimensional solutions. By substituting Eq.(4) into Eq.(3), we obtain 



{k-l)^-ik-l)b-—l 
dr m 



u{r) + 



(fu(r) f , m'\ du(r) m' , , ../ x , x ^, ^ , ^ 

^ ' ' 2b-\ — p^ H bu{r) - 2mV{r)u{r) + 2E{m - mo)u{r) 



r = 0. (5) 



dr"^ \ m J dr m 

Next we consider a series solutions for the above reduced radial wave Schrodinger equatim. 

III. THE SERIES SOLUTION WITH A CLASS OF STATIC POTENTIALS 

Equation (3) cannot be solved exactly except for some particular cases. Neverthless, we 
can get an approximate solution using the series expansion method [18,19]. In application, 
we consider here a group of central potentials belong to the following general form [23] 

V{r) = -Vir-" + V2r^ + V3, (6) 

where Vi and V2 are positive coupling constants whereas the constant V3 may be of ei- 
ther sign. Moreover, this group of potentials satisfies the boundary conditions stated in 
[20,21,22,23]. This class of generality for potentials (6) is used to produce the bound state 
energy spectra for quarkonium systems [20,21,22,23]. It comprises a well-known potential, 
e.g., the Cornell potential {we set a = (3 = 1, Vi = A = ^as, V2 = B, V3 = C) (cf. Ref.[21]). 
Now we try the following series expansions about origin: 

'nT'ir) = J2 ^1^'^'" = bo + J2 ^'^^''' ^0 = mo (7) 

m' °° 

together with the series expansion for u{r) 

00 
u{r)=Y,air\ ao 7^ (9) 

and substitute Eqs.(6)-(9) into Eq.(5), we obtain the following relation 

Y, [-{k - l)hair' - 2biair' + {k - l)iair'^^ + i{i - l)air'"^ 

i=0 

00 00 

- J2 K {iarr^^' + lay+' - bay+'+') + b^air'+^ + 2EJ2 Kair"-^'^' 

u=0 u=0 



oo 



u=0 u=0 u=0 

with 6^ = —2moE. On the other hand, we define 



uO-ir 



u+i+l 



(10) 



j+u=i 

Mi= J2 «i^-; ^>o (11) 

j,u=0 

M[=J2a,h',- 2>0 (12) 

Tr= JZ 3Ci,K- 2>1, To = 0. (13) 

Setting the coefficients of the power of r" to be zero, we obtain the following recurrence 
relation of the bound energy spectrum 

n{n + l)a„+i + (A; - l)(n + l)a„+i - {k - l)6a„ - 26na„ - IM'^ + 6M^_i - T„ 

+2^M„_i + 62a„„i + 2\/iM„+«_i - 2\/2M„_/3_i - 2\/3M„_i = 0, (14) 

with the final radial wave functions 

00 
RuAr) = r^^-^^'^e-'"- ^ a^ , a, ^ 0. (15) 

j=0 

We now present some special cases. 

A. Coulomb Potential 

For the solutions of the Coulomb problem in A^-dimensional space, we set Vi = Z, 
^ = V3 = 0, q; = 1, /3 = 0, then we get its recurrence relation from Eq.(14): 

n{n + l)a„+i + (/c — l)(n + l)an+i — {k — l)6a„ — 26na„ 

-/M; + 6M;;_i - T„ + 2EM„_i + 62a„_i + 2ZM„ = 0, (16) 

with the radial wave functions are given in Eq.(15), cf. Ref. [24]. This case was treated in 
Ref..[19]. 



B. Harmonic Oscillator Potential 

For the solutions of the harmonic oscillator problem in iV-dimensional space, we set 
V2 = uP' 1 Vi = V3 = 0, q; = 0, /3 = 2, then the recurrence relation (14) becomes 

n{n + l)a„+i + (A; — l)(n + l)a„+i — (k — l)6a„ — 2bnan 

-IM'„ + hM'^^^ - T„ + 2EMn-i + 6'a„_i - 2oo^Mr,^^ = 0, (17) 

with the radial wave functions are given in Eq.(15). This case was also treated in Ref.[19].. 

C. Confining Linear Potential 

For the solutions of the confining linear potential in A^-dimensional space, we set Vi = 
V3 = 0, V2 = B, a = 0, (3 = 1, then Eq.(14) becomes 

n{n + l)a„+i + {k — l){n + l)an+i — {k — l)6a„ — 2bnan — IM'^ 

+bM^_^ - T„ + 2EM„_i + 62a„„i - 25M„_2 = 0, (18) 

with the radial wave functions are still given in Eq.(15). 

D. Cornell Potential 

Now we investigate a confinement potential consisting of an attractive Coulomb term 
and a confining linear potential used for calculation of quarkonium (qq) bound-state masses 
[20,21]. With the set of parameters Vi = A = 4a,/3, V2 = B, V3 = C, a = (3 = 1, then the 
recurrence relation (14) becomes 

n{n + l)a„+i + (A; - l)(n + l)a„+i ~ {k - l)6a„ - 26na„ - IM'^ + hM'^_^ - Tn 

+2EM„_i + 6^a„_i + 2AM„ - 25M„_2 - 2CM„_i = 0, (19) 



with the radial wave functions are given in Eq.(15). 

Let us investigate the last case. The recurrence relation (19) implies 

Wq - 2Amo' 



ai 



k-1 



ao, 



(20) 



_ [{k + l)b + {I + 1% - 2Amo] 
^' 2~k ^' + 



[{k + 3)6 + (/ + 2% - 2Amo] 

^3 = -^ TTl 7T^ 0,2 + ^^ -^ — T -ai 



lb[ - bb'o - 


- 2Abi + 2Cmo 




2k 


[(/ + 1)6; - 


- bb'Q - 2Abi + 2Cmo 



3{k + l) 



3(A; + 1) 



(21) 



[/6'2 - bb[ + 2{C-E)bi- 2Ab2 + 2Bmo] 



(22) 



IV. THE ANALYTICAL SERIES SOLUTIONS FOR M = MqE-^^ 

We choose one simple example considered recently by Ref.[18]. Assuming a particle with 
an exponentially decaying position-dependent effective mass m{r) = moe~^^ , A > [18]. 
This form is taken on the base that the position-dependent mass must be convergent when 
r ^ oo. We try the series expansion 



m{r) = '^m^r", mo = 1, 

u=0 



(23) 



and 



m 



m 



-A. 



(24) 



Thus, for Cornell potential, Eq.(5) gives the following recursion relation 



[{k - l)n + n{n - 1)] a„ - [b{k - 1) + (26 - A)(n - 1) - /A] a„_i 



+6(6 - A)a„_2 + 2Em„_2a„_2 + 2Am„_ia„_i - 25m„_3a„_3 - 2Cm„_2a„_2 = 0, (25) 



which implies 



ttl 



b- 



l\ + 2Amo 
k-1 



Oo, 



(26) 



_ [{k + 1)6 - A(/ + 1) - 2Amo] [2Cmo + Xb - 2Ami] 

0-2 — ~i '^l "I TTi '^O) 



2k 



2k 



(27) 



[{k + 3)6 - (/ + 2)A - 2Amo] [Xb + 2Cmo - 2Ami] 

CI3 — ;77^ — —;-{ CI2 H ^77^ — —-{ Oi 



3(A; + 1) 



3(A;+1) 



[2Bmo + 2{C - E)mi - 2Am-^ 

H TTT, — —7^ CtO- 



(28) 



3(A; + 1) 

Finally, for another application of the model in case of a constant particle mass m = mo. 
We find the Coulomb's wave functions. in three-dimensional space, A^ = 3, easily through 
the relations (26)-(28). We present the following results [24] 

2Am() n 



Oi 



(n + / + 1) {21 + 2) 



Oo, 



(29) 



02 



n{n — 1) 



2Amo 
n + l + lj 2!(2/ + 3)(2/ + 2 



-Oo, 



(30) 



as 



2Amn 



n{n — l){n — 2) 



n + l + lj 3!(2/ + 4)(2/ + 3)(2/ + 2) 
and from which the wave functions is given by 



ttQ, 



2Amn 



'^ ' U, \n + l + lj {2l+l + i)\{n-i)\i\ 



(2/ + l)!n! 



ao. 



with a proper normalization 



'^0 ~ ^^n,l 



2 Am. ^'+^ 



n + / + l 



Am^oul 



-, n = 0,l,2,-- 



(31) 



(32) 



(33) 



(n + / + 1) V (n + 2/ + 1) 

On the other hand, considering a particle with an exponentially decaying position-dependent 
effective mass m,{r) = iriQe"^^ , the relations (26)-(28) read [24] 



Oi 



Am,Q 



i + i 



IX 



n + l + l 2moA 



flO, 



(34) 



Arrin 



a2 



{21 + 3) 



1 + 2 
n + l + l 



2Amo 



ai + 



A mi 

2(n + / + l) mo 



oo 



(35) 



as 



yl?72n 



3(/ + 2) 



1 + 3 {l + 2)\ 



n + l + 1 



2Amn 



- 1 



02 + 



A 



2(n + / + l) 



mi 
mo 



tti 



Ami ^^2 

2{n + l + 1) mo 

and from which the Coulomb's wave functions 



flO 



R^iir) = r'+^e T^^^^'' \ao + air + a2r^ + a^r^ + 



where a^, with i = 0,1, 2, 3, are given by means of Eqs.(33)-(36). 



V. CONCLUDING REMARKS 



(36) 



(37) 



We have carried out the series solutions of position-dependent mass A^-dimensional ra- 
dial Schrodinger equation for a general class of static potentials which are mostly used 
in the calculations of quarkonium bound-state masses. Some interesting results including 
the Coulomb, linear, harmonic, and Cornell potentials are also considered. In particu- 
lar, the Cornell case was investigated and the series solution for the A^— dimensional ra- 
dial Schrodinger equation with an exponentialy decaying position-dependent effective mass 
m,{r) = moc"'^'", are solved analytically. Contrary to most of the approximation methods 
which are valid when the parameters are small, our results obtained by this approximation 
method are also valid with large values of the parameters. These analytical calculations can 
be applied to some related fields of physics such as quarkonium systems in three dimensional 
space, i.e., A^ = 3. The work presented through the two examples discussed above would 
appear to be of interest in the ground that it offers an explicit formulas used to construct 
the wave functions for a central potential [20,21,22,23,24]. 
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